LGP-30 USERS' ORGANIZATION - POOL 
"SEARCH FOR ROOT" 
Program No. C2-37 

PROGRAM DESCRIPTION 

Program Title : "Search For Root" 

Author : . Arthur I. Larky and Jerry Rayna 

Installation : Lehigh University 

Purpcae : 

To solve for the real roots of equations of the form 

f (x) = 

by an iteration process. (The method of false position ), 

The program may be used either as a keyboard routine to print 
X and f (x) or as a subroutine within a program. 

Input: 

Two distinct initial guesses for each root must be stored in L Q 
+ 0051 and L + 00^2. 

A constant t stored in L + 0062 indicating the tolerance allowed 
in f (x). 

An accelerating factor may be required (optional). 
(See also Procedure). 

Output : 

When used as a subroutine only: X in the accumulator at a q 
determined by the subroutine used to calculate f(x), (See below). 

When used as a keyboard routine: 

Decimal printout of X and f (x). With transfer control depressed, 
the above is printed for each iteration. If the control button is up, 
printout will occur only for those values of f (x) which are less than 
€ in absolute value. 




R L Q + 0110 7 

U L r if X-i x and f- are to be inputed. 
o <- 

R L Q +■ 0110 "/ . „ v v J „ , J . 

U L + 0002 T 1* 2 an( * € are alreac fe r in memory 

The routine will exit with X in the accumulator. 



Hartree, "Numerical Analysis, " Oxford University Press, Second Edition 
(1958) Page 211;. _ 1 _ 
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As a keyboard routines 

Transfer to L 

Type in X,, X 2 , and 6 

into E + 0051, L Q + 0052, and L Q + 0062 using the 

proper data input identification words. 

Depress transfer control if every iteration is to be printed. 
Routine will continue iteration until X is printed twice after which 
an overflow stop will occur. ° 

Subroutines Required: 

1) A properly scaled subroutine capable of calculating f(x) for the 

X values most likely to be encountered must be in memory. The R-U calling 
sequence for this subroutine must be placed into L + 0003 - L + OOOlj. 
and L + 0009 - L Q + 0010. ° 

2) A data input subroutine with calling sequence in L and L + 0001 
(optional), 

3) A data output subroutine with calling sequences in L + oil8 - 0120 
and L Q + 0127 - 0129. (optional). 

h) A test quadratic E2x + 5x +1 - f(x)"7 is included in the coding 
and is called by the sequence R0113, U0101 iri locations 0003-1* and 0009-10. 
x l» x 2» and 6 must be entered @ q =» k as usual. 

Proc edure ; 

l) The iteration uses the formula: 

X i+1 = X i - 1 + x i - 1 " X i • f (xi-l) 



l_f (x ± ) _ f ( x . _ X )J / 

where = |f(x i )| + jf(x i _ 1 )| + 1 at 30 
and is repeated until j f (x) ) — ^ 



2) For use as a subroutine in a general program: The general pro- 
gram must provide X-^, X and € in the proper memory locations or read them 
from tape. The transfer control button may be depressed during program check- 
out to print X and f(x) for each iteration (except the final one). Routine 
exits with X in accumulator. 

3) For use as a keyboard routine: 

Change L Q + 0027 to U0028. Then X and f (x) will be printed 
(at q's determined by the data output calling sequence) when | f (x) I <£, ^ 
or for each iteration when the transfer control is depressed. 

h) An accelerating factor to aid in convergence may be included by 
changing L Q + 0155 to M0113. The accelerating factor (at q - 0) should be 
stored in L Q + 0113. (OPTIONAL). 

-2- 
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f>) Break Points 16 and 32 should be depressed. 

Progr amm ed Stops; 

Two break point 32 stops and one break point 16 stop-' are included. 
The breakpoint 32 stops follow the f(x) print and the X. + 1 calculation. 
Their purpose is to permit the iteration to be halted and new values for 

1» ^2> an< * € to be read in. This is done by skipping tbe execution of 
the next program step. To do so release the Break Point 32 button. When 
the computer stops, perform the following steps: 

a) depress One Operation button 

b) depress Manual Input button (on console) 

c) depress Start button 

d) return to Normal Operation and Start. 

The break point 16 stop enables the programmer to look at the size 
of the correction factor for the iteration. 

Limitations ; 

When used as a keyboard routine, the program will over flow stop 
after the second printout of the root. To search for additional roots, 
release the Break Point 32 button, and restart and when a break point 
32 stop occurs, return to data input as outlined under Programmed Stops 
above. * ~ 

It is likely that, if the initial guesses for X^ and X 2 are too far 
away from X q3 overflow stops may occur. In this case merely depress the 
start button as often as necessary. Regardless of the presence of an over- 
flow, some correction to X^ will be effected and ultimately values within 
scale will be obtained and the search continued. The function calculation 
subroutine must, of course, be able to handle X values reasonably near its 
root without overflow. The method causes looping at a point of inflection. 

Storage: 

12? locations, no temporary storage. 



-3- 
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f- 
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^2 aropp 
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■A ' 
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t 
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Hoot search Example : -nbirar 

2X 2 + 5X + 1 - f(X). 

1. With transfer contrux button 

.0000800' 

o»o4o85i ' y -0000005 ' .. ooo< jooc ' 
4+040862' 1' -ooocooc. " 



W** ^K**' - w rx n * ? 



uu: 



>>:■-. rro^raM No. C2- 



b? '.I. I^rky anrf 1. .Jayna 



Transfer tc .... 



The two guesses fur > 
The epailon. 



-2.28079589 

-2.280776^5 
-2.28077645 



i ■ K ■ 
.OU0O8O 

.00000c 

.000000 



Final Solution. 



O+040851V 0' -OOvOvua. ' xaa/Ka, ' 



•.2I85U294 
.. 21922728 
-.21922357 
..21922357 



.OO2807 
.000014 
.000000 
.000000 



Nev guesses for > 



F.nai Solution 



2. tfltffr transfer control button down: 

(Mote effect, itf any, of depression of transfer control out ton or. 
other subroutines ir, asej In particular, the lata input rout, nn . 



.0000800 s 

0+040651 ' 3' 
4+o4o862« 1* 



-OOCAAAJ*; 

- 0000000 : 



■ AK.Kf.*.f 



Transfer to ic. 

The two guesses; Lei X . 
The epailon. 



2_.999599_ 
-3.800000 

-7.000000 
-3. 144 578 

-2.6l4026 

-2.413957 

-2.306810 

-2.282337 
-2.280795 
-2.280776 
-2.280776 



f. \t.) 

^■9999 
010.8799 
064.0000 

005.0536 

002.7673 
OOG.5845 
000.1086 
J00.Q064 

tOO.OOOG 
000.0000 
0OG.000C 



Solution after first iteration. 
Solution after second iteration. 



Final Solution 



0+040651 ' 1 -0000001' -0000000 s 

001, c coo 
-000. 4444 

-ooc.047' 



0.000000 

-0.333333 

: p r 230769 

-'■■ ; - rOdS^ 

-O.S.V---S27 
-0.21.. ,.23 

-0.219223 



-000,0000 

ooo.ooco 

OOO.OCJOC 



Vew guasses for X. 

SoJ vfcio--. ifter first iteration. 



Final solution 



S" 
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"SEARCH ZZ-3. ROOT" 



PAGE OF 



2A/60 



TRACK 



PROGRAM INPUT CODES 



LOCATION 



jlOlo 











OjLfc 



FORM LP-IQ 
PRINTED iN U 



°±M, 



INSTRUCTION 



OPERATION 



ADDRESS 



111 



o i Ok 



1" 



Oil 



Lj 



-L_JL 



to il 



I 

I 



-4-- 



|2 



L1L 



16 



j9 



I I! 



|8 



2 |0 



2 I 



2 13 



2 |4 



2 15 



IL 



2 |S 



n. 



x 



R. U! 7lQ 



Hi 21 7 1 



Bl 0i0l5 



R| 0|1|1 



U| 01 HO 



H{ OlOlg 



U_| Oi 0|0 



B_jO_lOji 



R| Q| 111 



Ul Ol liO 



Hi 0| 0j6 



Bl Oi 0,5 



H|Q ,0,5 



U| 0| 0|2 



B| OlOig 



H| 0|0|5 



uj Q|Q|i 



_ B .j °l°l 6 



H 1 1 O il 



T j | nu 



P|0|0|2 



[ I I 



B j & Oil 



hi 0IQI5 



U| 0|0[1 



1 \ 0|°l 6 



T [0 | 1 10 



.j Ql 1 |i 



P| 01013 



+-»-*- 



CONTENTS 
OP ACTRESS 



NOTES 



* 



Xi 



S 



f(xi) 



V/T 



i ./ 



Xi + 1 



Data Input Routine 



Tor Xj,, X P , £ 



calc. f(xi) 



Z3 



f( xi+l) 



f(xi) 



.raw j'C ri-l ) 



lX3 Xi+1 . 



' I f (xi+l) 



j2 



E 



m 



N Xi 



H 



.£ 



ZEI 



S.A. 



Royal McBee Corporation x carriage return 

DATA PROCESSING DIV. — r 7 — 
if * 



IEI 



calc:. f(xi+l) 



new.Xi 



urv, f(yj) 



ch ange eign_ 



new Xi-1 



search over 



print Xi, f(xi) on but - 

ton 



f(xL-l)/0 



PORT CHESTER, NEW YORK 



/ = CONDITIONAL STOP CODE 
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JOB NO. 


PROGRAM NO. 

C2-37 


RROG*AM fSJffcARfcP 5XL j PROGRAM CHECKED BY: 

Gl Rayna FOOL Review 


°?h/60 


PROSLEW: "SEARCH FOR ROOT' 

j 


TRACK 


PROGRAM INPUT CODES 


a. 
O 

t/1 


LOCATION 


INSTRUCTION 


a. 
O 


CONTENTS 
CF ADDRESS 


NOTES 


OPERATION j ADDRESS 


1 
1 1 1 | 1 1 1 


1 


















X 








OlO |3 |2 


iFlO 1 OlO I 0|0|0 






$ @ 6 






1 13 |3 


1 1 U f._.I. ...I . . 1... 

I 1 |B|0 | 0|^|5 




f(xi) 








1 |3 |4 


1 J 1 1. J.. .-l....|._, 

I I |T|0 i 1|0,6 






corapl. 






1 13 |5 


1 ! 1 j. ...1-.„„J.-.„,I... . 

i 1 1 U|0 I 0,3|6 


>y\ __ 








1 P |6 


— i — i — i — , — i — i — i — i 

i 1 1 H|0 i 0|5|8 




f(xi) 1 








1 |3 V 


1 X -J 1 1 1 I - 

1 1 iB(0|O|£|6 




ftY-i-l) 








1 13 |8 


— i — i — ) — j .i — i -j. — 
I 1 |T|0 | 0|U|0 






COIMDl , 






1 13 |9 


! — ■* ' ' ! ' ' — ' 

! 1 1 UIOI 01 lil2 


'X- - 








1 |4 |0 


* — ■ — ' — r*- 4 — *■■** l ■ ■ 
1 1 I C|0| 1|0|5 






slear accumulator 






1 |4 |l 


1 " 

i 1 1 sjoi oi 516 




f(v\--\) 




1 J 1 ...-L__, 




1 H |2 


1 1 1 AIOI li 010 




1 © 30 








i I 4 |3 


*■'— J ' ■■■*■ ■■( "■ ■ ,J " 

1 1 1 AjOl 01 518 


liMiihtii 








1 |4 |4 


1 ! t H| i 1| 2 1 3 













1 |4 IS 


i i i u| on nn 








jiom ml m on I 9 


! 


1 |4 |6 


i i I I i > i 












1 |4 |7 


i i i 1 i i i 


■ ■■■ w 


X J 1 






1 |4 |8 


■■' ' I ' 












1 14 |9 


8| 0|0 | 0|0|0 | 0(0 




Cm 06) 


-1 @ 






1 15 10 


— ' — *——* — 1 — ' — "■■ — *■ 

1 1 i 1 1 1 1 












\ 15 |l 


— ' 1 1 — -{ »- — I — > — 

c. i i 1 i I .:/ 




X 


Y-i+T 






1 |S |2 


L i , i I i P 






Xi 






1 15 |3 


c . ! ■ , p 






Xi-1 






1 15 14 


'— ■ i ' ■ ' 

1 ! I 1 1 1 1 












1 15 |S 


c, i , i i i p 




X 


f(xi) 






1 |S |6 


■ i * 
£. i i 1 i i P 






ffxi-1) 






1 |5 |7 


i i i i i i i 








I I i | | i I 




1 15 |B 


C , ■ , , D 








f (*i) 1 






1 IS |9 


1 > 1 1 1 1 I 




X 








1 I* 10 


1 1 1 1 1 | 1 












1 16 II 


•C! , i , , P 






f (vi +1 ) 






1 |6 |2 


— 1 ' — 

U I 1 | 1 I J 






6 


i i i 1 1 I i 




| |6 |3 


1 

1 1 1 1 I 1 1 




X 




FORM LP- 12 
PRINTED IN U.S.A. 








Royal I 

DA 


4cBee Corporation* 

TA PROCESSING DIV. — fj " 


5 
/ 


V CARRIAGE K 
= CONDITIONA 


ETURN 
L STOP C 


:ooe 
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PROGRAM CHECKED BY: 

POOL Review 


DATE „ , . 

2A/60 


problem: "SEARCH FOR ROOT" 


TRACK 


PROGRAM INPUT CODES 


o 


LOCATION 


INSTRUCTION 


O 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS 


f <" 
i f I i l i i 


/ 














i 

1 1 ! I 1 ! | 




X 




i ~ — — 
< t ! I 1 | I 




0,1 io io 


i I 1 I 1 1 12 






1 ® ^0 


s — . — i — i — i 




1 |0 |l 


— i — i — i — j — i — i — i — 
1 1 IH |0I1 1 5|1 




x @ h 








1 10 |2 


I | IM 1 Oil 1 2th 




2 @ 2 




I ■ f I ! ! ! 




1 |0 |3 


I 1 1 Aj 010 i 312 




X 


$ ® 6 








1 1° I 4 


i | | Uj 1 1| 1| 1 












1 1° I 5 


G i , ■ , P 




(OOhO) 


dump 






j |0 |6 


_HEI i L ( 1 1..-1. . 

1 1 1 M] OIO 1 hi 9 




-1 & 








1 |0 |7 


! 1 lU 1 010 i 3l6 




x 








| |0 |8 


r — ' — ' — u i — ' — * ■"' 
1 ! i l l l 1 








i \ 1 | 1 1 1 




I 1° I 9 


i 1 1 B|0|0 | 5|2 




Xi 








1 |l |0 


1 1 1 UtQiO ( OJPJ 






"R" to here for search 






i 1' H 


■ — ' — ' — ' — p™- 1 - J — ,,rfr 

1 ! 1 MjOll 1 511 




X 


X @ h 








i |i 12 


1 1 1 AjO|1 1 612 




1 @ 10 








i |l !3 


i I | uj/pio i oioj 












1 II 14 


8 | 0| X | Fi 0|0 | 0|0 






-1 @ 






1 |l |5 


— i — i — i... _[...,. i... .i 
I IX 1 PM 1 6 )ilh 




[X 




C.R. 






t 1' I 6 


1 1 IB 1 010 1 512 




Xi 








1 P 1"' 


■—'- l ' 1 ■'■■■' ■ ' 

1 |X|F J 0|0 | 0|3 






dslav 






1 1' |8 


i IXiRI3I0i1|2 






"~\Data Output 






| l» I' 


I — ■ — ' — ' i r ' — L ■■■■*■ 

1 IXI Ul 310 I Oi 




X 




V 






I I 2 !° 


i 

i |Xi F10I0 ! Oil 






5 print xi 






l I 2 I 1 


.... ,,!,.. „l..... i. | 

i 1 lU I Oil 1 216 








^O'O I <">|n 1 n|n| }| 


t 


1 |2 12 


1 i I 1 i 1 1 












1 |2 |3 


C . . ! . . 3 




X 










1 (2 U 


hO | i OJ 0|0 | 0|0 






2 ® 2 






1 |2 I? 


C , , j . i P 






A*/V 






1 |2 |6 


t i |B|0|0|5|5 




fCvil 








1 |2 I? 


i ix m 1 310 i 11 2 




X 




^ Data Output 






1 |2 |8 


i 

l 1 X IU 1 310 i 01 






/ 






1 1 2 1 "? 


I 

I IX IF j 0| Oi Of 3 






print f(xi) 






1 |3 |0 


1 IX IF | 312 1 6|1 






Bk Pt 32 STOP 


> • i LI ,1 > 




_i ! 3 l« 


i i |U|0|0|3!3 




XI 




FORM LP-IO 
PRINTED IN U.5.A. 
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DA 
POR 
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TA, PROCESSING DIV. ~ V "' 
r CHESTER, NEW YORK 


/ 


E 


CARRIAGE t 
CONDITIONA 


ETURN 

I STOP CODE 
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? *$8 ff W i0 - 




PROGRAM CHECKED 3Y: 

POOL Review 


PATE 

2-U -60 


Pi!OBLEM: "SEARCH FOR ROOT" 


TRACK 


PROGRAM INPUT CODES 


fa- 

o 


LOCATION 


INSTRUCTION 


ex. 
O 

S; 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION | ADDRESS 


1 


/ 
















X 






0) 113 |2 


i t | U |0 1 tO 1 






read new X, £ 






1 |3 |3 


I I |B|0 |0,S| 5 




f(xi) 




> I 1 j I 1,1. 




1 13 |4 


I l 1 S i 0* 1 5 i 6 




f(xi-l) 








1 |J |5 


. ' " '" 'I'" ' "''■■■ 
i 1 IU| 0ili3 1 7 




X 








1 13 |6 


'" * " ■' — ' ■■"»■• • ■ " "' 

i 1 1 | I i 1 












1 P |7 


1 1 iDj 0,1,2| 3 




y 








I 13 |S 


i I |U| 0|li3 I 9 


' 










1 P |V 


I 1 iHj 01,2 i 5 


/ 


X 










1 |1 |0 


1 I ,u! 0, l,Uj 1 












1 1* |l 


I 1 .3(0,0,5,6 




f(xi-l) 








1 M 12 


, | |U|0 i i,U|U 












1 1* |3 


I 

I i ! | 1 1 1 


/ 


Aj 








1 H |4 


, | ,DjO,l|2,3 




X 








I 1* |5 


J 1 10(0 I 0131 1 


/ 


*r 




1 ! 1 ] ! 1 ! 




1 1* 16 


1 1 ! B I° I 0|5 f 3 




Xi-1 








1 h I? 


, | |Sj0| 0,5 1 2 


/ 


X 


JXj 




1 f I j | l 1 




1 |4 IS 


I 1 |HjO,1,5| 2 












... J J4.1'. 


1 1 |M|0|1|1|U 




-ion 








| |S |0 


I ' '" ■ 

l 1 i UjOl 0|2|6 












1 I 5 i' 


C i i | i i D 


1 ■ 


XI 


v ® 9 






.1 |5 J? 


' "'■ ■"'■ ■■?■■■ "'■'"" ■ 

, | , MfO| 0|3[1 




f(xl-l)/ 









1 IS (3 


" "1 " " ■'' 

l 1 I DJQI 11215 












1 |5 |4 


1 . XIF|1, 6|2 S 6 






Stop Bk. Pt. 16 


1 1 1 j 1 i 1 




! IS |5 


1 1 1 U|0| 1j5i6 




X3 




1 1 I I | 1 1 




1 |S |6 


i 1 lAjOl 01513 


..... , . . 

Ll Xi-1 








1 15 |7 


i I | H10I 0|5ll 


' 










1 IS |8 


I lX|F|3l2 j 310 






Stop Bk. Pt. 32 






_J._ I 5 i!_. 


i i 1 UjO|0 i 0|9 


/ 


X! 








1 . I« |0 


1 1 ! UjO|0 • 0.0 








^n i nin|n i n\n i ? 


1 


1 | b !' 


i 

i i i 1 i t 1 












1 |4 |2 


1 12 1 OjOl 010)0 






1 © 10 


III!!!. 




1 1* |3 


1 1 1 1 1 1 1 


/ 


y\ 




r-ORM LP-12 
PRINTED IN U.S.A. 
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BURN 

I STOP CODE 





